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ABSTRACT 

Brown dwarfs are important objects because they may provide a missing link between stars and planets, 
two populations that have dramatically different formation history. In this paper, we present the candidate 
binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered 
during 2004 - 201 1 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate 
events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG- 
399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-201 l-BLG-149, and MOA- 
201-BLG-278/OGLE-201 1-BLG-012N. Among them, we are able to confirm that the companions of the lenses 
of MOA-201 1-BLG-104/OGLE-201 l-BLG-0172 and MOA-201 l-BLG-149 are brown dwarfs by determining 
the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The 
measured mass of the brown dwarf companions are (0.02 ± 0.01) M Q and (0.019 ± 0.002) M Q for MOA-201 1- 
BLG-104/OGLE-201 l-BLG-0172 and MOA-201 l-BLG-149, respectively, and both companions are orbiting 
low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of 
lensing events with well covered light curves increases with new generation searches. 
Subject headings: gravitational lensing: micro - brown dwarfs - binaries: general 



1. INTRODUCTION much higher than planets. Studies of brown dwarfs are im- 

Brown dwarfs are sub-stellar objects that are too low in P ortant ]> e f use ™y Provide a missing link between 
mass to sustain hydrogen fusion reactions in their cores but stars and P lanets d Kulkarni| | 1997|) , two populations that have 
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dramatically different formation history. In addition, the 
Galaxy may be teeming with brown dwarfs although there 
is no evidence for a large population of brown dwarfs from 
current observations dGraff & Freesell99 6: Najita et al.ll2000t 
iTisserand et al.l f2007). Brown dwarfs had lon g been though t 
to exist based on theoretical considerations (lKumarl[l969l) . 
However, these objects are intrinsically faint and thus the 
first confirmed brown dwarf was n ot discovered until 1995 
( Reb olo et al]fl99l lNakaiimal[l995l) . 

There have been 2 major methods of detecting brown 
dwarfs. The first is direct imaging by using ground (e.g., 
Schuh 2003) and space-borne (e.g., Mainzer 2011) infrared 
instruments. As in the case of extrasolar planets, brown 
dwarfs can also be indirectly discovered by detecting wob- 
bles in the motion of the companion star or the flux decrease 
of the companion star caused by the brown dwarf occulting 
the companion star. 

Microlensing is also an effective method to detect brown 
dwarfs. Gravitational microlensing refers to the astronomical 
phenomenon wherein the brightness of a star is magnified by 
the bending of light caused by the gravity of an intervening 
object (lens) between the background star (source) and an ob- 
server. Since the phenomenon occurs regardless of the lens 
brightness, microlensing was proposed to detect dark compo- 
nents of the Ga laxy such as b lack holes, neutron stars, and 
brown dwarfs (Paczvriski 1986). 

Although effective, the application of the microlensing 
method in brown dwarf detections has been limited. The main 
reason for this is the difficulty in distinguishing brown dwarf 
events from those produced by main sequence stars. For gen- 
eral lensing events in which a single mass causes the brighten- 
ing of a background star, the magnification of the lensed star 
flux depends on the projected lens-source separation, u, by 



u 2 + 2 



M ( M 2 +4)V2' 



(1) 



where the lens-source separation is normalized by the radius 
of the Einstein ring, 6>e (Einstein radius). The lens-source sep- 
aration is related to the lensing parameters of the time scale 
for the source to cross the Einstein radius, f E (Einstein time 
scale), the time of the closest lens-source approach, to, and the 
lens-source separation at the moment of the closest approach, 
uo (impact parameter), by 
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Among these lensing parameters, only the Einstein time scale 
provides information about the physical parameters of the 
lens. However, it results from the combination of the lens 
mass, distance, and transverse speed of the relative lens- 
source motion and thus it is difficult to uniquely determine 
the mass of the lensing object. Not knowing the mass, then, 
it is difficult to single out brown dwarf events from those pro- 
duced by stars. In principle, it is possible to measure the lens 
mass by additionally measuring the Einstein radius and the 
lens parallax, which are respectively related to the physical 
parameters of the lens by 



9 E = (KMTr Iei ) 
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and 
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where 7r rel = AU(D L ' -Dg 1 ), K = 4G/(c 2 AU), M is the mass 
of the lens, and Z\ and £>s are the distances to the lens and 
source star, respectively. The Einstein radius is measured by 
analyzing the distortion of the l ensing light curve affected by 
the finite size of the source star (iNemiroff & W ickramasinghe 
ll994HWitt&Maolll994HGouldlll994h . while the lens paral 
lax is measured by analyzing the deviation in a lensing light 
curve caused by the deviation of the relative lens-source mo- 
tion from a rectilinear one due to the change of the observer's 
positi on induced by the orbital motion of the Earth around the 
Sun ((Gould 1992). Unfortunately, the chance of simultane- 
ously measuring 6*e and tte is very low for a single-lens event. 
The finite source effect is important only for very high mag- 
nification events, for whic h the lens-source s eparation is com- 
parable to the source size dChoi et alj|2012l) . but these events 
are very rare. The parallax effect is important only for events 
with time scales that are a significant fraction of the orbital 
period of the Earth, i.e. 1 yr. However, such long time-scale 
events are unlikely to be produced by low mass brown dwarfs 
because f E oc M l l 2 . As a result, there exist only two cases for 
which the le ns of a single lensing event was iden tified as a 
brown dwarf (ISmifh et alJl2"003l:lGould et al.ll2009l) [ T l. 

In binary lensing events, in which the lens is composed of 
two masses, on the other hand, the chance to identify the lens 
as a brown dwarf is relatively high. There are several rea- 
sons for this. First, for binary lensing events, it is possible 
to routinely measure the mass ratio, q, of the lens compo- 
nents. Then, brown dwarf candidates can be sorted out based 
on the first order criterion of small mass ratios. Second, the 
majority of binary lensing events involves causti c crossings 
dSchneider & Weil3ll986l:lMao & Paczvnskil 19911) . The caus- 
tics represent the positions on the source plane at which the 
lensing magnification of a point source becomes infinite. The 
caustic-crossing part of a lensing light curve varies depending 
on the source size. Analysis of the light curve affected by this 
finite source effect allows one to measure the Einstein radius, 
which allows for better constraint on the lens mass. Third, the 
chance to measure the lens parallax is high, too. This is be- 
cause the duration of an event corresponds to the total mass of 
the binary lens system not the low mass component and thus 
the average time scale of a binary lens event is longer than 
that of a single lens event. 

Despite these advantages, there exists only one event 
(OGLE-2008-BLG-5 10) for which the companion of a binary 
lens i s known to be a brown dwarf candidate (Boz za et alj 
l20T2h . The main reason for the rarity of brown dwarf events 
is the difficulty in analyzing binary lensing light curves due to 
the complexity of the light curves. For binary lensing events 
in which the mass ratio between the lens components is very 
low as in the case of a binary composed of a star and a planet, 
the signature of the low mass companion appears as a short 
term deviation on the top of an otherwise smooth and symmet- 
ric single lensing light curve of the primary. As a result, iden- 
tification and analysis of planetary events are relatively sim- 
ple, resulting in more than 20 identified microlensing plan- 
ets. However, the perturbative nature of the companion breaks 
down and the light curve becomes very complex for binaries 
with mass ratios greater than q ~ 0.1, which corresponds to 
the mass ratio of a typical binary lens composed of a brown 

1 Among these two cases, the microlensing brown dwarf discovered by 
Gould et al. 1 2009) is exception that provides the rule because the lens mass 
could be determined not by measuring the lens parallax induced by the 
Earth's orbital motion but from the "terrestrial parallax". 
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dwarf and a star. This introduces difficulties in the immediate 
identification of a brown dwarf event. A firm identification 
of a brown dwarf event requires detailed modeling of lensing 
light curves. Unfortunately, modeling light curves of binary 
lensing events is a difficult task due to the complexity of \ 2 
surface in the parameter space combined with the heavy com- 
putation required for analysis. As a result, detailed modeling 
of binary events was conducted for a small subset of events, 
resulting in a small number of identified brown dwarfs. In re- 
cent years, however, there has been great progress in the anal- 
ysis of binary events with the development of analysis meth- 
ods based on advanced logic combined with improved com- 
puting power. This progress made the analysis not only more 
precise but also faster. As a result, nearly all binary lensing 
events are routinely modeled in real time wit h the progress of 
events in current microlensing experiments ( iRvu et aLll2010b 
iBozza et alJ2012h . 

In this paper, we present brown dwarf candidates found 
from the analyses of binary lensing events detected from mi- 
crolensing experiments conducted from 2004 to 201 1. In §2, 
we describe the criteria of choosing brown dwarf candidates, 
observation of events, and data reduction. In §3, we describe 
the procedure of modeling. In §4, we present the results of 
the analysis. In §5, we draw our conclusions. 

2. SAMPLE AND OBSERVATIONS 

To search for brown dwarf candidates, we investigate bi- 
nary lensing events discovered during 2004 - 201 1 microlens- 
ing observation seasons. Considering that typical Galac- 
tic microlensing ev ents are produced by low mass stars 
dHan & Go uld 2003), the lower mass companion of an event 
with a mass ratio q < 0.2 is likely to be a brown dwarf can- 
didate. Therefore, we sort out brown dwarf candidate events 
with the criterion of q < 0.2. For this, we conduct modeling 
of all binary lensing events with relatively good light curve 
coverage combined with obvious binary lensing features such 
as the spikes caused by crossings over the fold or approaches 
to the cusp of a caustic. The good coverage criterion is 
needed to secure the mass ratio measurement. The criterion 
of the characteristic features is needed because it is known 
that binary lensing events with weak signals can often be in- 
terpreted as other types of ano malies such as binary sources 
(Jaroszvhski & Skowron 2008). For binary lensing events 
detected before 20 08, we also refer to the previous anal- 
yses c onduct ed by [j aroszvhs ki et all (120061) . ISkowron et al.l 
(120071) . and llaroszyriski et al.1 (l2010h . From the 2004 - 
2011 search, we found 7 candidate events, including OGLE- 
2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG- 
006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA- 
201 1-BLG-104/OGLE-201 l-BLG-0172, MOA-201 1-BLG- 
149, and MOA-201 -BLG-278/OGLE-201 1-BLG-012N. 

In Table [T] we list the candidate events that are ana- 
lyzed in this work along with the equatorial and Galactic 
coordinates of the lensed stars. The name of each lens- 
ing event is formed by taking the survey group who dis- 
covered the event (i.e. OGLE or MOA), followed by the 
year of discovery (e.g. 2004), the direction of the observa- 
tion field ("BLG" designating the Galactic bulge field), and 
a sequential number assigned to the event. If an event is 
discovered independently by two different groups, they are 
named separately. Survey observations were conducted by 
the two groups: the Optical Gravitational Lensing Experiment 
(OGLE: Udalski 2003) and the Microlensing Observation in 
Astrophysics (MOA: Bond et al. 2001, Sumi et al. 2003). 
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FIG. 1.— Light curve of OGLE-2004-BLG-035. The notation in the paren- 
theses after the label of the observatory denotes the pass band of observation. 
The insets show the enlargement of the light curve during caustic crossings 
and approaches. 
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FIG. 2.— Light curve of OGLE-2004-BLG-039. Notations are same as in 
FigQ] 
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FIG. 3.— Light curve of OGLE-2007-BLG-006. Notations are same as in 
FigfT] The notation for the observed passband "N" denotes no filter was used. 

Among 7 events, 5 were additionall y observed by follow-up 
observation groups in cluding uFU N (lGouldll2006[). PL ANET 
dBeaulieu et all 12006b. RoboNet (fTsapras et all (20091) . and 
MiNDSTEp (iDominik etail 12010b . In Table 12 we list the 
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TABLE 1 
Coordinates of Events 



event 



OGLE-2004-BLG-035 

OGLE-2004-BLG-039 

OGLE-2007-BLG-006 

OGLE-2007-BLG-399/MOA-2007-BLG-334 

MOA-201 1-BLG-104/OGLE-201 l-BLG-0172 

MOA-2011-BLG-149 

MOA-201 1-BLG-278/OGLE-201 1-BLG-012N 



TABLE 2 
Telescopes 



event telescopes 

OGLE-2004-BLG-035 OGLE, 1.3 m Warsaw, Las Campanas, Chile 
OGLE-2004-BLG-039 OGLE, 1.3 m Warsaw, Las Campanas, Chile 
OGLE-2007-BLG-006 OGLE, 1.3 m Warsaw, Las Campanas, Chile 
/LtFUN, 1.3 m SMARTS, CTIO, Chile 
/^FUN, 0.4 m Auckland, New Zealand 

/JUN, 0.4 m FCO, New Zealand 

OGLE-2007-BLG-399 OGLE, 1.3 m Warsaw, Las Campanas, Chile 
/MOA-2007-BLG-334 MOA, 2.0 m Mt. John, New Zealand 
AtFUN, 1.3 m SMARTS, CTIO, Chile 
PLANET, 1.0 m SAAO, South Africa 

MOA-201 l-BLG-104 OGLE, 1.3 m Warsaw, Las Campanas, Chile 

/OGLE-201 l-BLG-0172 MOA, 2.0 m Mt. John, New Zealand 
AtFUN, 1.3 m SMARTS, CTIO, Chile 
/^FUN, 0.4 m Auckland, New Zealand 
/uFUN, 0.4 m FCO, New Zealand 
AtFUN, 0.5 m Wise, Israel 
/uFUN, 0.3 m PEST, Australia 
RoboNet, 2.0 m FTN, Hawaii 

RoboNet, 2.0 m FTS, Australia 

MOA-201 l-BLG-149 MOA, 2.0 m Mt. John, New Zealand 

OGLE, 1.3 m Warsaw, Las Campanas, Chile 
AtFUN, 1.3 m SMARTS, CTIO, Chile 
/xFUN, 0.4 m Auckland, New Zealand 
RoboNet, 2.0 m FTS, Australia 

MiNDSTEp, 1.54 m Danish, La Silla, Chile 

MOA-201 l-BLG-278 MOA, 2.0 m Mt. John, New Zealand 
/OGLE-201 1-BLG-012N OGLE, 1.3 m Warsaw, Las Campanas, Chile 
/iFUN, 1.3 m SMARTS, CTIO, Chile 
/^FUN, 0.4 m Auckland, New Zealand 
/xFUN, 0.4 m FCO, New Zealand 
/^FUN, 0.4 m Kumeu, New Zealand 
/^FUN, 1.0 m Lemmon, Arizona 
RoboNet, 2.0 m FTS, Australia 
MiNDSTEp, 1.54 m Danish, La Silla, Chile 

NOTE. — CTIO: Cerro Tololo Inter-American Observatory; FCO: Farm 
Cove Observatory; SAAO: South Africa Astronomical Observatory; PEST: 
Perth Extrasolar Survey Telescope; FTN: Faulkes North; FTS: Faulkes 
South 

survey and follow-up groups who participated in the obser- 
vation of the individual events along with the telescopes and 
their locations. 

Photometric reductions were carried out using photome- 
try codes developed by the individual groups. The OGLE 
and MOA data were redu ced by the photometry c odes de- 
veloped by lUdalskil (120031) and iBond etalJ (120011) , respec- 
tively, whic h are based on the D ifference Image Analy- 
sis method (lAlard & Lup ton 1998). The fiFUN data were 
processed using a pi peline based on the DoPHOT software 
dSchechter et all [19931) . For PLANET and MiNDSTEp da ta, 
a pipeline based on the pySIS software ( Albr ow et alj2 009) is 
used. For the red uction of RoboNet data, the DanDIA pipeline 
(Bramich 2008) is used. Photometry errors of the individual 
data sets were rescaled so that x 2 P er degree of freedom be- 
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comes unity for the data set of each observatory, where % 2 
is estimated based on the best fit model. We eliminate data 
points with large errors and those lying beyond 3a from the 
best fit model to minimize their effect on modeling. In Figures 
1 - 7, we present the light curves of the events. 

3. MODELING 

In the analyses of each lensing light curve, we search for 
the solution of lensing parameters that yield the best fit to the 
observed light curve. The basic structure of a binary lens- 
ing light curve is described by 7 parameters. The first three 
parameters to, uq, and fg are identical to those of a single- 
lensing event and they specify the source motion with re- 
spect to the lens. The next 3 parameters describe the binary 
lens system, including the projected separation between lens 
components in units of the Einstein radius, s, the mass ra- 
tio between the lens components, q, and the angle between 
source trajectory and the binary axis, a. The final param- 
eter is the source radius, 6L, normalized to the Einstein ra- 
dius, p+ = 0+/0E (normalized source radius). The normalized 
source radius is needed to account for dilution of the lensing 
magnification caused by the finite size of the source. The fi- 
nite source effect is important when the source is located at 
a position where the magnification varies very rapidly such 
as near the caustic and thus different parts of the source star 
are magnified by different amounts. In addition to these basic 
binary lensing parameters, it is often needed to consider the 
parallax effect for a fraction of events with long time scales. 
To describe the lens parallax effect, it is required to include 
the two additional parameters, he.n and he.e, which represent 
the two components of the lens parallax vector tte projected 
on the sky along the north and east equatorial coordinates, re- 
spectively. The size of the lens parallax vector corresponds to 
the ratio of the Earth's orbit to the Einstein radius projected on 
the observer observer's plane, i.e. 7r E = AU/|>eZ)s/(Ds - ,Dl)], 
where r E = £>l#e is the physical size of the Einstein radius. 
See Equation (4). 

In the analysis, we search for the best fit solution in two 
steps. In the first step, grid search is conducted for a subset 
of parameters while o ther parameters ar e optimized by using 
a downhill approach (Dongetal. 2006). By inspecting the 
X 2 distribution in the grid parameter space, we identify lo- 
cal minima. In the second step, we investigate the individual 
local minima and further refine the solutions by allowing all 
parameters to vary. Once the solutions corresponding to the 
local minima are found, we choose the best fit solution by 
comparing the x 2 values of the individual local minima. In 
this process, we also identify possible degenerate solutions 
resulting from various causes. We set s, q, and a as grid pa- 
rameters because these parameters are related to light curve 
features in a complex way such that a small change in the 
values of the parameters can lead to dramatic changes in the 
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FIG. 4.— Light curve of OGLE-2007-BLG-399/MOA-2007-BLG-334. 
Notations are same as in FigfT] 
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FIG. 5.— Light curve of MOA-201 1-BLG-104/OGLE-201 l-BLG-0172. 
Notations are same as in Fig|3] 
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FIG. 6. — Light curve of MOA-201 l-BLG-149. Notations are same as in 
FigLI] 

resulting light curve. For the x 2 minimization in the downhill 
approach, we use the Markov Chain Monte Carlo (MCMC) 
method. Once the solution of the parameters is found, we es- 
timate the uncertainties of the individual parameters based on 
the distributions of the parameters obtained from the MCMC 
chain of solutions. 
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FIG. 7.— Light curve of MOA-201 1-BLG-278/OGLE-201 1-BLG-012N. 
Notations are same as in Fig|3] 



We use the inv erse ray-shooting technique 
( Schn eider & Weifil 119861: iKayser et all 1 19861 IWambsganssI 
119971) to compute magnifications of events affected by 
the finite source effect. In this numerical method, rays 
are uniformly shot from the image plane, bent by the lens 
equation, and land on the source plane. Then, the finite 
magnification is computed as the ratio of the number density 
of rays on the source plane to that on the image plane. The 
lens equation for a binary lens, which describes the relation 
between the positions of a ray on the image and source 
planes, is represented by 



m 2 



Z-ZLA Z~Zl,2 



(5) 



where m\ = \/(\ + q) and /«2 = qm\ are the mass fractions 
of the binary lens components, £ = £ + ir\, z = x + iy, and 
Zl = *l + iyh represent the complex notations of the source, 
image, and lens locations, respectively, and z denotes the 
complex conjugate of z- The ray-shooting method requires 
heavy computation. We thus minimize computation time 
by applying the ray-shooting method only when the source 
crosses or approaches very close to caustics and utilizing 
the semi-analytic hexadec apole approximation dGouldll2008l : 
iPeicha & Hevro vskv 2009) in the adjacent area. 

In computing finite magnifications, we consider the varia- 
tion of the magnification caused by the limb-darkening of the 
source star surface for events with well-resolved features in- 
volving caustics. In our model, the surface brightness profile 
is described by 



Sx = 



Fx 



l-Tx 1 



- COS if) 



(6) 



where Fx is the linear limb-darkening coefficient, F\ is the 
source star flux, and if) is the angle between the normal to 
the source star's surface and the line of sight toward the star. 
We set the limb-darkening coefficients corresponding to the 
source type that are determined based on the color and mag- 
nitude of the source. In Table[3] we present the adopted limb- 
darkening coefficients, corresponding source types, and the 
measured de-reddened color along with the assumed values 
of the effective temperature, T e a, the surface turbulence ve- 
locity, v tur b, and the surface gravity, logg. 

For all 7 events, we conduct modeling based on the 7 ba- 
sic binary lensing parameters. We carry out additional mod- 
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FIG. 8. — Geometry of the binary lens systems for the light curves presented in Fig. 1 - 7. In each panel, the two dots (marked by Mi and Mi) represent the 
locations of the binary lens components, where the bigger (Mi) and smaller (M2) dots represent the higher and lower mass binary components, respectively. 
The figure composed of concave curves is the caustic. The line with an arrow represents the source trajectory. For three events of MOA-201 1-BLG- 104/OGLE- 
2011-BLG-0172, MOA-201 l-BLG-149 and MOA-201 1-BLG-278/OGLE-201 1-BLG-012N, we present the geometry of two solutions caused by the close/wide 
degeneracy. To better show the caustic-involved perturbations region for each of the three events with degenerate close and wide solutions, we present the 
enlargement of the perturbation region (enclosed by a box) in a separate panel on the right side of each main panel. 



eling considering the parallax effect for events with Einstein 
time scales fg > 30 days. When the parallax effect is con- 
sidered, we separately investigate solutions with impact pa- 
rameters uq > and mq < 0. This is because the lensing light 
curves resulting from the source trajectories with +uq and -mo 
are different due to the asymmetry of the trajectory caused by 
the parallax effect, while the light curves are identical when 
the trajectories are straight lines due to the symmetry of the 
magnification pattern with respect to the binary axis. 

4. RESULTS 

In Table HI we present the solutions of the parameters de- 
termined from modeling for the individual events. For events 
where the additional parallax modeling is conducted, we also 
list the parallax solution. In addition, we present the solu- 
tions of local minima if they exist. The model light curves 
resulting from the best fit solutions of the individual events 
are presented on the top of the observed data points in Figures 
1-7. In Figure [8j we also present the geometry of the lens 
system corresponding to the best fit solution of each event. 
Based on the best fit solutions, it is found that the mass ra- 
tios of the analyzed events are in the range 0.09 < q < 0.20, 
indicating that the lower mass companions of the lenses are 



strong brown dwarf candidates. 

For 3 out of the total 7 events, it is found that there exist lo- 
cal minima other than the best fit solution. The events with 
degenerate solutions include MOA-201 1-BLG- 104/OGLE- 
2011-BLG-0172, MOA-201 l-BLG-149, and MOA-201 1- 
BLG-278/OGLE-2011-BLG-012N. For all of these events, 
the degeneracy is caused by the symmetry of the lens equa- 
tions between the close {s < 1) and wide (s > 1) bina- 
ries which is kn own as the "close/wide binary degeneracy" 
dDominikl IT9991) . For MOA-201 1-BLG- 104/OGLE-201 1- 
BLG-0172 and MOA-201 l-BLG-149, the degeneracy is very 
severe with A\ 2 = 2.5 and 3.1, respectively. For MOA-201 1- 
BLG-278/OGLE-201 1-BLG-012N, the degeneracy is less se- 
vere and the wide binary solution is preferred over the close 
binary solution with A\ 2 = 259.4. We note that the mass ra- 
tios of the pair of close and wide binary solutions are simi- 
lar to each other although the binary separations are very dif- 
ferent. Therefore, the degeneracy does not affect the brown 
dwarf candidacy of the lens. 

For three events, it is possible to measure the lens mass by 
simultaneously measuring the lens parallax and the Einstein 
radius. With these values, the lens mass and distance are de- 
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TABLE 3 

Limb-darkening coefficients and source information 



quantity 


MOA-2011-BLG-104 
/OGLE-2011-BLG-0172 


MOA-2011-BLG-149 


MOA-2011-BLG-278 
/OGLE-2011-BLG-012N 


IV 


0.51 


0.73 


0.59 


r« 


0.43 


0.64 


0.51 


Ti 


0.36 


0.53 


0.43 


(V-I)o 


0.62 


1.64 


0.85 


source type 


FV 


KV 


GV 


Tiff (K) 


6650 


4410 


5790 


Vmrb (kms ') 


2 


2 


2 


logg (cms -2 ) 


4.5 


4.5 


4.5 



TABLE 4 
Model Parameters 



event 


X 2 /dof 


to (HJD') 


lift 


tr (days) 


s 


q 


Cl 


P* (10~ 3 ) 




7TJ7 £ 


OGLE-2004-BLG-035 


270.8/240 


3085.825 


0.1502 


66.81 


1.021 


0.122 


2.129 


2.36 


- 


- 






±0.090 


±0.0011 


±0.29 


±0.002 


±0.002 


±0.014 


±0.16 








246.6/240 


3085.651 


0.1475 


65.00 


1.024 


0.130 


2.113 


2.24 


-0.02 


0.07 






±0.167 


±0.0017 


±0.68 


±0.003 


±0.003 


±0.019 


±0.27 


±0.06 


±0.02 


OGLE-2004-BLG-039 


393.0/386 


3081.422 


0.1009 


38.95 


1.163 


0.093 


6.188 


6.64 










±0.054 


±0.0031 


±0.67 


±0.005 


±0.004 


±0.008 


±0.39 






OGLE-2007-BLG-006 


325.6/587 


4171.109 


0.0043 


60.92 


0.801 


0.110 


6.219 


1.70 










±0.161 


±0.0003 


±3.53 


±0.014 


±0.009 


±0.004 


±0.15 






OGLE-2007-BLG-399 


1289.8/999 


4305.037 


0.1906 


23.39 


1.446 


0.168 


1.855 








/MOA-2007-BLG-334 




±0.129 


±0.0244 


±2.23 


±0.047 


±0.016 


±0.017 








MOA-2011-BLG-104 


1726.9/3210 


5670.304 


0.0478 


38.86 


0.379 


0.089 


1.638 


1.50 






/OGLE-2011-BLG-0172 




±0.005 


±0.0008 


±0.48 


±0.004 


±0.003 


±0.003 


±0.12 








1724.6/3210 


5670.310 


0.0478 


40.06 


3.213 


0.118 


1.639 


1.36 










±0.005 


±0.0007 


±0.48 


±0.048 


±0.005 


±0.003 


±0.13 








1712.4/3208 


5670.286 


-0.0481 


38.47 


0.379 


0.090 


-1.628 


1.54 


-0.19 


0.10 






±0.006 


±0.0009 


±0.49 


±0.005 


±0.003 


±0.003 


±0.13 


±0.30 


±0.03 




1709.9/3208 


5670.291 


0.0494 


39.29 


3.195 


0.121 


1.628 


1.36 


-0.33 


0.08 






±0.006 


±0.0008 


±0.52 


±0.048 


±0.005 


±0.003 


±0.13 


±0.21 


±0.03 


MOA-2011-BLG-149 


5031.1/4918 


5700.402 


0.0081 


148.60 


0.256 


0.145 


2.450 


0.55 










±0.008 


±0.0003 


±4.97 


±0.003 


±0.005 


±0.004 


±0.02 








5036.1/4918 


5700.447 


0.0080 


150.60 


5.077 


0.229 


2.464 


0.55 










±0.008 


±0.0003 


±5.74 


±0.082 


±0.013 


±0.004 


±0.02 








4889.0/4916 


5700.399 


0.0067 


179.73 


0.240 


0.134 


2.446 


0.44 


-0.74 


-0.24 






±0.008 


±0.0003 


±8.50 


±0.004 


±0.007 


±0.004 


±0.02 


±0.04 


±0.02 




4892.1/4916 


5700.441 


0.0064 


183.91 


5.210 


0.190 


2.460 


0.44 


-0.69 


-0.20 






±0.008 


±0.0002 


±5.66 


±0.088 


±0.012 


±0.003 


±0.02 


±0.04 


±0.02 


MOA-2011-BLG-278 


5820.7/5558 


5744.741 


0.0128 


18.96 


0.265 


0.130 


0.532 


2.77 






/OGLE-2011-BLG-012N 




±0.001 


±0.0002 


±0.26 


±0.001 


±0.003 


±0.001 


±0.04 








5561.3/5558 


5744.731 


0.0140 


17.12 


4.597 


0.200 


0.514 


3.12 










±0.001 


±0.0002 


±0.21 


±0.027 


±0.005 


±0.001 


±0.04 







NOTE. — HJD'=HJD-24500()0. 

termined by 

Be 

M=— , (7) 

K7T E 

and 



7rEf?E + 7TS ' 

respectively. Here its = AU/Ds represents the parallax of the 
source star. 

For MOA-2011-BLG-104/OGLE-2011-BLG-0172, the 
value of the lens parallax measured from modeling is 
7r E = (7r E W 2 + 7r E £ 2 ) 1 ' /2 = 0.34 ±0.21. The Einstein radius is 
measured based on the normalized source radius determined 
from modeling combined with the information of the angular 
source size by 9% = #*/ p+. The angular source size is assessed 
from the de-reddened color and brightness of the source star, 
which are estimated based on the source brightness measured 
in V and / pass bands and calibrated by the relative position 



. MOA-201.1-BL(3.'-104--. 
- ■ /OGLE-2011-.BL&.--0JP 
giant clump. -:.: . ' ■.' 
I .source . ■ Vfv" '..^ 

"-'--V :■■ -^Wi- ' 



mi 



-1 -0.5 
(7-7) 




FIG. 9. — Color-magnitude diagram of stars in the field of MOA-2011- 
BLG-104/OGLE-2011-BLG-0172 and MOA-201 l-BLG-149 with the loca- 
tions of the lensed stars and giant clumps marked as blue and red dots, re- 
spectively. 



of the source on the color-magnitude diagram with respect to 
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a reference (Yoo et al. 2004). We use the centroid of bulge gi- 
ant clump as the reference because its de-redd ened brightness 
of In, p = 14.45 at the Galactocentric distance dNataf & Gould! 
20120 and the color of (V-/) 0)C = 1.06 dBensbv et al.ll20Tlh 
are known. We measure the de-reddened color and brightness 
of the source star to be (V-I,I) , S = (0.70,18.15). Figure 
[9] shows the locations of the source star and the centroid 
of the giant clump in the color-magnitude diagram that 
is constructed based on the CTIO data. After the initial 
photometric characterization of the source, we learned that 
spectroscopic observation is conducted when the source was 
highly magnified (Bensby 2012, private communication). 
From this spectroscopic observation, the de-reddened color 
of the source star is estimated as (V—T)o,s = 0.62, which is 
adopted for our analysis. Once the de-reddened V — I color of 
the source star is measured, it is translated i nto V — K color 
by usi ng the V — I versus V — K relations of iBessell & Brett! 
( 1988) and then the source angular radius is estimated by us- 
ing the r elation between the V— K color and the angular radius 
given bv lKervella et alj d2004l) . The measured angular source 
radius and the Einstein radius are 9+ = (0.68 ±0.06) /ias and 
9e = (0.54 ± 0.07) mas, respectively. With the measured 
Einstein radius, the mass and distance to the lens are esti- 
mated as M = (0.20 ± 0.12) M© and D L = (3.29 ± 1 .20) kpc, 
respectively. Then, the corresponding masses of the primary 
and the companion are Mi =M/(l+q) = (0.18 ±0.11) M© 
and M 2 = qM\ = (0.02 ±0.01) M©, respectively. Therefore, 
the binary lens is composed of a low mass M dwarf and a 
brown dwarf. 

For MOA-2011-BLG-149, the lens parallax is also mea- 
sured with 7Te = 0.78 ± 0.04. For this event, however, multi- 
band (V and /) observations taken at a single observatory 
(CTIO) covered only near the high magnification peak of the 
light curve. As a result, the source brightness in both bands 
cannot be estimated by the usual method based on modeling. 
Fortunately, there exists no significant blending for this event 
and thus the color measured at the peak can be approximated 
as the source color. From the source brightness and magnifi- 
cation at the peak, we also estimate the baseline source bright- 
ness. The source location estimated in this way is plotted 
on the color-magnitude diagram presented in Figure [9] Once 
the source location is determined, the de-reddened color and 
brightness are determined by the usual way of using the clump 
centroid. The measured values are (V— I,T)o, s = (1.64,20.8), 
which correspond to those of a Galactic bulge K-type dwarf 
with an angular radius of 9+ = (0.46 ± 0.04) /ias. Then, the 
estimated Einstein radius is 9% = (1.04 ±0.10) mas and the 
resulting lens mass and distance are M = (0.16 ±0.02) M© 
and £>l = (1.07 ±0.10) kpc, respectively. The masses of the 
individual lens components are Mi = (0.14 ±0.02) M© and 
M 2 = (0.019 ±0.002) M©, respectively. Therefore, the lens 
system is also composed of a low mass M dwarf and a brown 
dwarf. 

For another parallax event OGLE-2004-BLG-035, no 
multi-band observation is conducted. However, the source 
star of the event is very bright, implying that the source flux 
is not much affected by blended light. Under this approxima- 
tion, we measure the source color in the OGLE catalog. Based 
on this approximate color combined with the giant clump lo- 

2 For other Galactic bulge fields with Galactic coordinates (l,b), the I- 
band brightness is given by /o, c = 14.45—5 log(cos/ + sin/cos</>/.vm<?>), where 
<f> ~ 45° is the orientation angle of the Galactic bulge with respect to the line 
of sight. 



cation, we estimate the de-reddened color and brightness of 
the source star of (V— 7,/)o, s = (1.19, 15.83), which correspond 
to the values of a low luminosity giant star with a source ra- 
dius of 9+ = (3.72 ± 0.32) /ias. Considering that there exists 
a non-zero blended flux fraction Fb/F tot ~ 20%, we decrease 
the source radius by ~ (Fb/F tol )/2 = 10%. The obtained Ein- 
stein radius is 9e = (1.40 ±0.21) mas. Then, the mass of the 
lens system is M = (2.48 ± 0.87) M© and thus the mass of the 
companion is M 2 = (0.29 ± 0.10) M©. The companion mass 
is greater than the upper mass limit of brown dwarfs of 0.08 
M©. 

5. CONCLUSION 

We searched for candidate binaries with brown dwarf com- 
panions by analyzing binary microlensing events discovered 
during the 2004 - 2011 observing seasons. Under the low 
mass ratio criterion of q < 0.2, we found 7 candidate events, 
including OGLE-2004-BLG-035, OGLE-2004-BLG-039, 
OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007- 
BLG-3 34, MOA-20 11-BLG-l 04/OGLE-20 1 1 -BLG-0 172, 
MOA-2011-BLG-149, and MOA-20 1-BLG-278/OGLE- 
2011-BLG-012N. Among them, we confirmed that the 
companions of MOA-20 1 1 -BLG- 1 04/OGLE-20 1 1 -BLG- 
0172 and MOA-201 l-BLG-149 were brown dwarfs by 
measuring the lens masses. 

The number of microlensing brown dwarfs is expected to 
increase with the improvement of lensing surveys and analy- 
sis technique. With the adoption of a new wide field camera, 
the number of lensing events detected by the OGLE survey is 
increased by a factor more than 2. In addition, a new survey 
based on a network of 3 telescopes located at three different 
locations of the Southern Hemisphere (Korea Microlensing 
Telescope Network) is planned to operate from the 2014 sea- 
son. Along with the improvement in the observational side, 
the analysis technique has greatly improved for the last several 
years, enabling prompt and precise analysis of a large number 
of anomalous events. With this improvement, the number of 
microlensing brown dwarfs will rapidly increase, making mi- 
crolensing one of the major methods in the discovery of brown 
dwarfs. 
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